NASA Grant NAG 1 


541 


/ /£? aj<?c ty 

//j -/r - c r /€L 

/// 



FINAL REPORT 
v/ 


VIBRATIONS AND STRUCTUREBORNE NOISE IN SPACE STATION 


R. Vaicaitis, C.S. Lyrintzis and D. A. Bofilios 


Department of Civil Engineering and Engineering Mechanics 
Institute of Flight Structures 
Columbia University 
New York, NY 10027 


October 9, 1987 


(NASA-CR-181381) VIBRATIONS AND N87-29590 

STRUCTUREBORNE NOISE IN SPACE STATION Fiudl 
Report (Columbia Univ. ) 152 p Avail: NTIS 

HC A03/MF A 01 CS CL 22B Unclas 

G3/18 0102711 



TABLE OF CONTENTS 


SECTION PAGE 

1.0 SUMMARY 2 

2.0 INTRODUCTION 3 

3.0 ANALYTICAL MODELS 7 

3.1 Noise Generation Inside a Cylindrical 

Enclosure 7 

3.2 Noise Generation Inside A Rectangular 

Enclosure 10 

3.3 Dynamic Response of A Cylindrical Shell 13 

3.4 Response of a Double Wall Cylindrical Shell 16 

3.5 Response of Double Wall Composite Shells 18 

3.6 Response of a Stiffened Shell 19 

3.7 Vibration of Double Wall Circular Plates 20 

3.8 Response of Stiffened and Interconnected 

Structures 26 

3.8.1 Response of a Stiffened Panel to 

Random Loads 26 

3.8.2 Response of a Stiffened Beam 31 

3.8.3 Coupled Stiffened Panel-Beam 

Structure 32 

4.0 NUMERICAL RESULTS 42 

4.1 Noise Transmission into Cylindrical 

Enclosure 43 

4.1.1 Single Wall Aluminum Shell 43 

4.1.2 Double Wall Aluminum Shell 44 

4.1.3 Double Wall Composite Shell 46 


l 



50 


4.1.4 Discretely Stiffened Cylindrical 

Shell 

4.1.5 End Plates 54 

4.1.6 Total Interior Noise 56 

4.2 Structural Response and Noise Transmission 

of Rectangular Enclosures 57 

4.2.1 Response of Stiffened Panels to 

Point Forces and Point Couples 59 

4.2.2 Response of Stiffened Panel to 

Calculated Bending Moment Inputs 60 

4.2.3 Response of a Coupled Skin-Stringer- 

Beam Structure 61 

4.2.4 Noise Transmission 62 

4.2.5 Noise Transmission - A Parametric 

Study 65 

5.0 CONCLUSIONS 68 

6.0 REFERENCES 71 

FIGURES 76 


1 1 



FOREWORD 


This study was performed for NASA Langley Research Center, 
ANRD by the Department of Civil Engineering and Engineering 
Mechanics, Institute of Flight Structures, Columbia University 
under Grant NAG1-541. The NASA Langley technical monitor was Dr. 
A.C. Powell. 


1 



1.0 SUMMARY 


Analytical models have been developed to predict vibrations 
and structureborne noise generation of cylindrical and 
rectangular acoustic enclosures. These models are then used to 
determine structural vibration levels and interior noise to 
random point input forces. The guidelines developed in this 
study could provide preliminary information on acoustical and 
vibrational environments in space station habitability modules 
under orbital operations. The structural models include single 
wall monocoque shell, double wall shell, stiffened orthotropic 
shell, descretely stiffened flat panels, and a coupled system 
composed of a cantilever beam structure and a stiffened 


sidewall. Aluminum and fiber reinforced composite materials are 
considered for single and double wall shells. The end caps of 
the cylindrical enclosures are modeled either as single or double 
wall circular plates. Sound generation in the interior space is 


calculated by coupling the structui 


aCOu StlC 


field in the enclosure. Modal methods and transfer matrix 


techniques are used to obtain structural vibrations. Parametric 


studies are performed to determine the sensitivity of interior 
noise environment to changes in input, geometric and structural 


conditions . 
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2.0 INTRODUCTION 


Structural vibrations and noise could play a significant 
role on the quality of physical environment in space station for 
prolonged orbital operations. Even though it is not expected 
that noise will be a major risk factor for physical health, it 
could be a significant source of annoyance, speech interference 
and fatique to individuals working under these conditions [1]. 
Furthermore, excessive vibrations and noise could have an adverse 
effect on delicate scientific experiments and manufacturing 
processes that are being proposed under zero gravity conditions 
[2], To satisfy the various proposed design objectives, the 
manned space station will need to perform many mission support 
functions such as pressurized and unpressurized laboratories, 
base for attached payloads, communications, command and control 
support, assembly, deployment and construction, maintenance, 
servicing, life support systems, etc. Mechanical vibrations 
resulting from power supply units, life support systems, 
electrical equipment, control thruster action, etc., could induce 
unwanted vibrations and noise. The information that is presently 
available on vibro-acoust ic environment for space station 
operations is very limited. The experience from Skylab [3-7] and 
other sources [1,2] could serve as a starting point of 
identifying some of the potential vibration and noise problems in 
the space station. The structureborne noise induced by 
mechanical vibrations could be a significant contributor to the 
noise environment inside the habitability modules. A detailed 
account of structureborne noise related work was presented in a 
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review article [8] . In general, the generation and transmission 
of structureborne noise are not well understood and fundamental 
theoretical and experimental work is needed. 

This report presents an analytical study on vibration 
response and noise generation in cylindrical shells and 
rectangular enclosures due to mechanical random point loads. The 
geometry of these structures are taken to be representative of 
pressurized habitability modules of the space station design. 

The end caps of the cylindrical shell are modeled as circular 
plates. Single wall, double wall and stiffened (frames and 
stringers) shells are taken as different options of the 
structural model. For the double wall construction the exterior 
shell could serve as a radiation and/or thermal shield while the 
interior shell is the main load carrying structure. The space 
between the two thin shells is assumed to be filled with soft 
thermal insulating materials. For the low frequency range 
considered in this study, the stiffened shell is represented by 
- an equivalent orthotropic shell wherein the effect of frames and 
stringers is smeared into an equivalent skin [9,10], The 
rectangular acoustic enclosure model is used to approximate the 
interior space of the habitability modules and to study 
structureborne noise generation by vibration of small panels, 
partitions, stiffened panels and complex structural geometries 
involving stiffened beams and discretely stiffened panels. The 
space station artist's concept and a typical habitability module 
are shown in Figs. 1 and 2. The structural models considered in 
the present study are presented in Figs. 3-7. The loads are 
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random point forces which can be acting at any arbitrary location 
on the shell/ end plates, panels, stiffeners or the sub-structure 
as shown in Figs. 6 and 7. 

The shell skins are modeled according to thin shell theory 
[11-16] and the end plates using the theory of thin circular 
plates [17-24] . For the case of fiber reinforced composite 
materials, the equations of motion are developed for the cases 
where a shell is a composite buildup of laminae, which consists 
of fibers imbedded in a supporting matrix [25-30] . The solution 
for shell and end plate vibrations are developed utilizing a 
modal approach [31,32]. Similarly, the acoustic wave equation 
for the generated interior acoustic pressure due to shell and/or 
end plate vibrations is solved by the Galerkin-1 ike procedure. 
Hard walls and absorbent boundary conditions at the interior 
surface are considered. To determine vibration response and 
noise transmission for the interconnected structures shown in 
Fig. 7, a modified transfer matrix procedure was developed [33— 
35]. In this approach, arbitrary point loads and/or distributed 
loads can be acting on the structure (beam) and the main 
structure (stiffened sidewall). Due to the complexity of the 
load transfer paths of these built-up structures, it is not easy 
to construct models which couple the various sub-components into 
a single dynamic system. The transfer matrix method proved to be 
relatively straightforward to apply to these stiffened 
structures. However, procedures based on transfer matrices 
suffer some drawbacks in practice. Because of the successive 
matrix multiplication required in this approach, ill-conditioned 
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systems are produced which require high precision to obtain 
meaningful results. To circumvent these difficulties, several 
different methods have been used. These include the modified 
transfer matrix method [36.], delta matrix formulation [37,38], 
precalculation of products using symbolic manipulation algorithms 
[39], double and quadruple precision, and special purpose 
routines for manipulating large numbers. 
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3.0 ANALYTICAL MODELS 


3.1 Noise Generation Inside a Cylindrical Enclosure 

Consider a closed cylindrical enclosure with interior volume 
2 

V = n R L as shown in Figs. 3-6. The noise inside this 
enclosure is generated by the vibrations of the shell structure 
and/or the vibrations of the circular end plates. Assume the 
shell and the end plate motions are independent. The solutions 
for shell and plate vibrations due to random point loads are 
presented in Sec. 3.3 - 3.5. The acoustic pressure inside the 
cylindrical enclosure can be obtained from 

p = Pl + p 2 (1) 

where p-^ and P 2 are the acoustic pressures generated by the 
shell and end plate motions. The pressure p satisfies the 
linearized acoustic wave equation 


2 •• 2 

V p - pp = p/c ( 2 ) 


in which (3 and c are the acoustic damping and speed of sound 
inside the enclosure, a dot indicates time derivative, and 

222 22222 
V = a /dr + ( 1/r ) d/dr + ( 1/r ) a /59 + d /bx (3) 
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The boundary conditions to be satisfied are 

9pi/9r = -pw - ( p/Z ) i at r=R (4) 


9pi/9x = 0 at x = 0,L (5) 

9p 2 /9x = pw L + ( p/Z L ) p 2 at x = 0 (6) 

9p 2 /9x = -pw R - (p/Z R ) p 2 at x = L (7) 

9p 2 /9r = 0 at r = R (8) 


where p is air density , w, w L , w R are the displacements in the 
normal direction (positive outwards) of the shell, left end plate 
and right end plate, Z, Z L , and Z R are the absorbent wall 
impedances at the surfaces of the shell, left end plate and right 
end plate, respectively. In the present study, it was assumed 
that all interior walls are treated uniformly with insulating 
materials for which [40] 


Z = Z L = Z R = - pc {(1 + 0.0571 (2ttR 1 /pw) 0,754 ) (9) 

+ i ( 0.087 (2irR 1 /pw) 0 * 732 ) } 


where R-^ is the flow resistivity of porous acoustic material, 
a) is radial frequency and i = ^ -1. The acoustic damping in 
the interior is modeled as 
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( 10 ) 



where Jj is the Bessel function of the first kind of order j, 

= a jk //R where is the kth root of the equation 

dJj/dr = 0 . Substituting Eqs. 11 and 12 into Eqs. 2-8 and 
using orthogonality condition of acoustic modes gives a set of 
ordinary differential equations which can be solved for the 
generalized coordinates P^j and A detailed procedure of 

this approach is given in Refs. 31, 41,42. Then, the spectral 
density of the acoustic pressure p can be determined from 
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S (x,r,0,u) = S (x,r,9,o>) + S (x,r f 9,w) 
P Pi P 2 


( 15 ) 


where S and S are spectral densities corresponding to 
Pi P2 

acoustic pressures p^ and P 2 « The sound pressure inside the 
enclosure is calculated from 


2 

SPL(x, r, 9, go) = 10 log {S (x,r, 9 ,m) Am/pq) (16) 

P 


_ 9 

where pg is the reference pressure (p^ = 2.9 x 10 psi 
2 

= 20 nN/m ) • 


3.2 Noise Generation Inside a Rectangular Enclosure 

Consider a rectangular acoustic space occupying a volume V 
abd as shown in Fig. 7. Noise is generated in the acoustic 
enclosure through vibrations of the flexible portions of the 
sidewalls, partitions, or individual small panels which can be 
located at any arbitrary position on the structure. The 
perturbation pressure p within the enclosure satisfies Eq. 2 
where now 

2 2 2 2 2 2 2 
V = 9 / 9x + 9 /9y + 9 /9z (16) 

The types of boundary conditions to be satisfied by Eg. 2 depend 
on the interior surface conditions of the walls. These could 
range from acoustically hard walls to those of highly absorbent 
walls which are treated with acoustic insulation materials. At 
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acoustically rigid boundaries, 


9p/9n = 0 (17) 

where n is outward normal to the boundary. Equation 2 together 
with Eq. 17 can be used to calculate modes and modal frequencies 
in the enclosure. For a surface treated with acoustic absorbent 
materials, the boundary conditions can be represented' by 

9p/9n = -pp/Z (18) 

where Z is the point impedance defined in Eq. 9. For a flexible 
elastic surface treated with absorbing materials, the boundary 
conditions to be satisfied are 

9p/'9n = -pw - pjb/Z (19) 

where w is the normal displacement of the vibrating sidewall, 
partition, or a localized panel. Expressing the acoustic 
pressure in terms of orthogonal modes corresponding to hard walls 
at x=0,b and y=0,a we write 

CO CD 

p(x,y,z,w) = l l A. . (z,w)Y. . (x,y) (20) 

i=0 j = 0 13 130 

in which A^j are the acoustic modal coefficients and j are the 
acoustic hard wall modes 
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Y..,(x,y,z) = cos(i-rt/t>) cos(jny/a) cos(kuz/d) (21) 

1 j K 

The acoustic modal frequencies are 

2 2 2 l /o 

to.., = c[(i-n:/b) + ( j u/a ) + (kn/d) r 2 (22) 

13K 

The solution for the perturbation perssure p can be obtained by 
taking Fourier transformation of Eqs. 2,18 and 19, substituting 
Eq. 20 into these equations and using the orthogonality condition 
of the acoustic modes. Then, the sound pressure levels inside 
the enclosure can be calculated using Eq. 16. Noise generated in 
the interior by other vibrating surfaces located at x=0,b and/or 
z=0,d can be estimated by a similar procedure. Then, the total 
acoustic pressure p inside the enclosure can be obtained from 

P = Pi + P2 + • • • P n (23) 

where p^ ,P 2 , . . . ,p n are the acoustic pressure contributions from 
the n vibrating flexible parts surrounding the enclosure. If 
these pressure contributions are assumed to be independent, the 
spectral density of the total acoustic pressure p can be 
calculated from 


Sp ( x , y , z , «) 


S 

Pi 


( x ,y , z , to) 


+ S 


P 2 


( x ,y , z , to) 


(24) 


+ 


• « • 


S ( x , y , z , to) 
p n 
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A detailed treatment on noise generation inside rectangular 
enclosures can be found in Refs. 34,35,38 and 43. 

The solutions for interior acoustic pressure, p, and the 
sound pressure levels are functions of the structural vibrations 
at the boundaries of the acoustic space. Next, we obtain 
solutions for vibration response of cylindrical shells, circular 
plates and discretely stiffened rectangular panels. 


3.3 Dynamic Response of a Cylindrical Shell 

Figure 3 shows a cylindrical shell exposed to external 
and/or internal random point forces. The input forces are taken 
as stationary and Gaussian random processes which can be located 
at any arbitrary position on the shell. A Dirac delta function 
is used to define the location of point load. The external p e , 
and the internal p 1 , random loads are expressed in terms of two 
point forces and F 2 as [26,31,32] 

p e ( x , 0 , t ) = (1/At A® ) (F® (t) 6(x-xt ) 6(9 - ef ) (25) 
+ F 2 (t) 5 ( x-x 2 ) 6(9 - 0 2 ) > 

p 1 (x , 9 , t ) = (I/A 1 A 2 ) {F 1 ( t ) 6( x - x\) 6(9 - 9i) (26) 

+ F 2 ( t ) 5 ( x - xj) 6(9 - 02)} 

where the superscripts e and i denote the external and the 
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1 


internal loads, 6 is the Dirac delta function, and for the 
cylindrical shell shown in Fig. 3, A® = 1, A® = R + h, = 
and A 2 = R . The point loads are assumed to be independent and 
each characterized by a spectral density. For a space station 
operation, the point loads could be generated through various 
mounts and attachments by vibrations of mechanical and electrical 
equipment, thruster action and other mechanical impacts. 

Using the Donnel-Mushtari approximations for thin shells 
[44,45], the equations of motion for the outward normal 
displacement w can be written as 

8 2 4 4 4 .. 4^1 

DV w + (Eh/R ) 5 w/6x + p h V w + yti = V {p e - p 1 > (27) 

s 


where 


3 2 


D = 

an /1 2 ( 1 - v ) 



(28) 

4 

4 4 2 

CM 

CM 

4 

4 4 

V 

= 5 /ax + (2/R ) 

a /ax ae 

+ (1/R ) 

a /ae (29) 

8 

8 8 2 

8 6 2 

4 

8 4 4 

V 

= a /ax + ( 4 /r ) 

a /ax ae 

+ ( 6/R ) 

a /ax ae (30) 


6 

8 2 6 

8 

8 8 


+ (4/r ) a 

/ax ae + 

(1/R ) 5 

/ae 


in which E, p , v and y are modulus of elasticity, material 
s 

density, Poisson's ratio and viscous damping coefficient, 
respectively. The general solution to Eq. 27 can be obtained in 
terms of simply supported shell modes 
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(31) 


w(x,e,t) = l I 


m=l n=0 


A ( t ) X* (x,9) 
mn mn 


where X^ n (x,0) = s i n ( ra ^x/L ) cosn 0 . The input loads p e and p 1 
are also expanded in terms of the shell modes. Substitution of 
Eq. 31 into Eq. 27 and use of the orthogonality principle, gives 
a set of ordinary differential equations in A jnn . Solving these 
equations and utilizing the theory of random processes [46] 
spectral densities of the displacement response are 


S (x,9,o) - l £ J 1 
m=l n=0 r=l s=0 


S («) X X 
mnrs mn rs 


(32) 


where S mnrg are the cross spectral densities of the generalized 
coordinates A mn . 


s = H ( o) 

mnrs mn 


* 

H ( w) 
rs 


r — G 

Lb 

mnrs 


+ S 1 • ] 
mnrs 


[ is ) 


The frequency response function H mn and the generalized random 


forces S are 

mnrs 


H 


mn 


p h [go - iwy/p h - go 
K s mn ~ f / K s 


(34) 


= (S^ (a>) X (x®, e®) X (xl, 01 ) 

mnrs F i mn 1 1 rs 1 1 


.e „e 


(35) 


+ S F 2 (t0) X mn (x 2' 0 2) x j- s ( x 2 1 e%)}/R 


where S and S® are the spectral densities of the random 


Fl 


F 2 
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point forces F^(t) and F 2 (t). Similar expression can be written 

for the generalized random forces S 1 acting on the inner side 

mnrs 

of the shell. The natural frequencies of a cylindrical 

shell can be obtained from 

2 2 2 4 

= {(D/p s h) [(mrc/L) + (n/R) ] 

+ (E/p R ) (m-rc/L) }/ [(mi/L) (36) 

fa 

2 2 

+ (n/R) ] 

The solution for shell vibrations presented in this section can 
be used in Eq. 4 to obtain noise generated inside a cylindrical 
enclosure . 


3.4 Response of a Double Wall Cylindrical Shell 

A double wall shell construction shown in Fig. 4 can be used 

to represent a design where the interior shell is the pressurized 

module and the exterior shell is used as a radiation and thermal 

shield. Following Refs. 31,32,41,42 and using Eq. 27, the shell 

vibrations can be modeled by two coupled partial differential 

equations for normal deflections w r and w_ as 

L X 


D E V E W E + (E E h E /R E )d ‘* W E /dX ‘ l + P E 7 E^E + V F t : k * <’ * P -” T > 


E E"E ' ' E l "s ' "E "I 


( 37 ) 


+ c E w E + ( 1/3 )n> s w E + ( 1/6 )m g w I } = 


16 


V E p e (x, 6, t) 



( 38 ) 


D I 7 I W I + ( E i h i/ R j ) S^Wj/dx 4 + + v x i k s {w i” W E^ 

+ CjWj + (l/3)m s w 1 + (1/6 )m g w E } = - Vj P 1 (x,0 / t) 

where m s and k g are mass per unit area and core stiffness. The 
subscripts E and I denote the external and the internal shells, 
respectively. The viscoelastic core separating the two shells is 
taken to be relatively soft, so that bending and shearing 
stresses can be neglected, and subsequently the soft core is 
characterized by a uniaxial constitutive law. Such a model 
allows in-phase (flexural) and out-of-phase (dilatational ) 
motions of the double wall systems. The vibrations of the inner 
shell generate noise inside the cylindrical enclosure. However, 
the motions of the two shells are coupled through the core. 

Thus, vibrations of the outer shell could induce motions of the 
inner shell and consequently generate noise in the interior. 

The equations of motion of the double wall shells are solved 
by modal expansion methods. The solution of Eqs. 37 and 38 is 
expressed in terms of simply supported shell modes 

00 00 

w (x,0,t) = l l A*L X® (x,9) (39) 

h i L ’ mn mn 

m=l n=0 

00 ® t c 

w T (x,e,t) = l l AX (x,9) (40) 

I ' ^ mn mn 

m=l n=0 

E I 

where A and A are the generalized coordinates of external 
mn mn 

• s 

and internal shells, and X are the shell modes. The solution 

mn 

procedure for shell displacements w and w T are similar to the 

h 1 

procedure presented in Sec. 3.3. The details of the response 
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analysis of double wall shells to random point loads are given in 
Refs. 32 and 42. 


3.5 Response of Double Wall Composite Shells 

The design of space structures is impacted by the 
interactions of functional requirements, such as strength, 
stiffness, weight, reliability, etc. To accomodate many of these 
requirements, new design concepts for lower weight, extended 
service life and improved integrity are needed. It has been 
demonstrated that composite materials could give weight and 
structural integrity advantages over many commonly used materials 
[47-50] . However, the low-weight composites might not provide 
any advantage with respect to less response or reduced noise 
transmission. Past studies have demonstrated that sandwich 
constructions might be an efficient way of dissipating 
vibrational energy [51-54]. Thus, to satisfy the required 
vibroacous t ic environment, designs utilizing composite materials 
might need to be modified by including the double-wall sandwich 
concepts . 

The sandwich shell system is composed of two simply 
supported cylindrical shells and a soft viscoelastic core as 
shown in Fig. 4. Each shell is a composite buildup of laminae, 
which consists of fibers imbedded in a supporting matrix. The 
laminae can be oriented in any arbitrary direction. The fibers 
are basically the load carriers. The equations of motion are 
derived using assumptions similar to those given in Refs. 26,44 


18 



and 45. The viscoelastic core separating the two composite 
shells is taken to be relatively soft, so that bending and 
shearing stresses in the core can be neglected. The natural 
frequencies and vibration response are obtained for simply 
supported cases by modal solutions and a Galerkin-like 
procedure. The details of the theoretical formulation and 
analysis of the double wall composite shells are given in Refs. 
32 and 42. 


3.6 Response of a Stiffened Shell 

The design configuration of the habitability modules is 
expected to be a discretely stiffened cylindrical shell as shown 
in Fig. 5. Analytical formulations and response calculations 
have been performed for the cases where the stiffened shell shown 
in Fig. 5 is represented by an equivalent orthotropic shell. In 
this case the effect of rings and stiffeners is smeared into an 
equivalent skin. Then, the natural frequencies can be calculated 
by the procedures presented in Ref. 9 and 10. For the 
application of low frequency vibrations and noise transmission, 
such a model might be adequate to evaluate the noise criteria 
inside the habitability modules. 
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3.7 Vibration of Double Wall Circular Plates 

Consider the two circular plates shown in Figs. 4 and 6 are 
simply supported at the edges. The governing equations of motion 
of the two plates, coupled- through a linear soft core, can be 
written as [17-24,42] 

4 . •• p p- 

D t V w t + c T w T + m T w T + k s (w T - w g ) + (l/3)m s w T 

(41) 

+ ( 1/ 6 ) m g w g = - p T (r,0,t) 

VS + C B*B + ffi B W B + k S (W B “ V + (1/3)m s*B 

(42) 

+ ( 1/ 6 ) m g w j = p B (r,0,t) 

where 


D T,B E T,B h T,B' /12( 1 ~ V T,B J 


(43) 


m T,B ~ P T,B h T,B 


(44) 


P - n P h P 

m - p n 
s 's s 


(45) 


,4 , ,_*£ + i » + 1,-4) ( 

5r 2 r 5r r 2 


JL + 1 + _i j £) 

6r 2 r £r r 2 26 2 


(46) 


The subscripts T,B denote the top and bottom plates and s 

T B 

denotes the core. The loads p (r,9,t) and p (r,9,t) are the 
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random excitations applied to the top and bottom plates. In 
obtaining Eqs. 41 and 42 it was assumed that the mass of the core 
follows an apportioned linear distribution. The w are the 
normal displacements of the midsurfaces of the top (exterior) and 
the bottom (interior) circular plates. The boundary conditions 
to be satisfied are 


W T f B (r ' 9) = ° 


at r = R p 


(47) 


2 

5 w S w 

f + V (1 T ' B , 1 T,B, n . 

.B 1 5r 2 T,B ( r Sr ~2 “775 — 1 = 0 at r = 


. 2 

5 w r 


2 2 
r 56 


(48) 


The solution to Eqs. 41 and 42 can be expressed in terms of 
normal modes 


w T ( r f 0 , t ) 
w B (r,9,t) 


l 


s=0 


oo 


l 

q=l 


I 


s=0 


00 


I 

q-i 


T^sq* 


B A sq (t) 



( r , 0) 
(r, 0) 


(49) 

(50) 


where T A sq 
(exterior ) 

X sq (r ' 9) 


and g A sq are the 9 enera li ze d coordinates of top 
and bottom (interior) circular plates, and 
are the circular plate modes given by 


x gq(r,9) = R sq (r) cos 


R sq < r > ■ J s (k sq r) - 77^7 

s v sq ; 


(51) 

(52) 
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in which J s and I s are Bessel functions and modified Bessel 
functions of the first kind respectively, and X® is the qth 
root of the frequency equation 


J s+i a) , 'sn' 1 ' 2x < 53 > 

" J s (k) ' T~JTJ — 


Results given in Eq. 53 are obtained by substituting Eqs. 51 and 

52 into Eqs. 47 and 48 and using relationships which relate the 

derivatives of Bessel functions to higher order functions [24 and 

P 4 * 

65] . In Eq. 53 \ = kR , k = ■ and consequently 


X 


s 

sq 


k R 1 
sq 


(54) 


k 


4 

sq 


T,B 



m T , B^°T , B 


(55) 


Substituting Eqs. 49 and 50 into Eqs. 41 and 42 and using 
the orthogonality principle, gives a set of coupled differential 
equations in T A gq and B A sq . Taking the Fourier transform of 
these equations it can be shown that 


T A sq (u) = H S q (w) { B A sg (u,) (k s + d/6)m u> 2 ) 


(55) 


+ T P sq (u) /5 S qJ /m T 


B A sq (t0) = H sq (aj ) { T A sq (^) (k g 


+ ( 1/ 6 )m g <j ) 


( 56 ) 


+ B P sq {u)/ 3 sq }/m B 
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(57) 


H sg (w) 1/7 , B^sq “ “ Y T,B /m T,B + ~ WC T,B /m T,B + k s /m T,B^ 


Y T,B m T,B + (1 /3)m s 


(58) 


R P 2 it 


t R^o (aj) = + / / p T,B (r,e,w) (r,9)rdrd9 

1 f D SQ g g Sq 


(59) 


R P 2it , *°sq if 3 * 0 

Ogn = f / {X P (r,9)} 2 rdrd9 = { 

Sq 0 0 sq 2nQ Qq if s = 0 


_ (R ) 


Px2 


(Jc 2 (0 + (1 - 


•) j«<Ol 


(60) 


sq l “s '''sq' ' '* ,.s .2 ' “s' sq 

s sq 

- ^ l^^sg'Vl 1 ^ 1 + J S <X sq> I s + l (X lq ) > (61) 

X s ( \ „ ) 

o rr CH 

*-> m w "n 

+ (R P ) 2 J s (X sq } r ( , + s 2 . 2. s _ '2s ,, 

2 2 s 1(1 + s 2 /I s U sq ) T s vX sq y} 

s sq ' sq ' (62) 

rp q W k ^rp n / m T l D 

T,B sq sq T,B T,B 


Furthermore, a ( ) ' indicates differentiation with respect to the 

spatial variable r and a bar indicates transformed quantity. 

The excitations applied to the top and/or bottom circular 
plates are assumed to be uniform random pressure or random point 
loads as shown in Fig. 6 for which the spectral densities are 
specified. In the case of uniform pressure input the generalized 
random forces reduce to 
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T,B B sq (u)) = { 


« -T ; 6 i i 

2 n^ 0 q p 


s = 0 
s * 0 


(63) 


where 


, = J? (j a 5 ) . J o u Qq ) T # x s 

•oa k 0g (j i'V I()(x . h'V 1 


(64) 


_T g 

and p ' (u) is the Fourier transform for spatially uniform 

T B 

pressure input p ' (r,0,t) . 

The random loads acting on the top and bottom plates are 

T B T B 

expressed in terms of two point loads F^' and ? i' as 


p T (r,e,t) = (l/AjAj){Fj(t) 6(0-0 1 ) 6(r-r x ) 
+ F T ( t ) 6(0-62) 6( r-r^) } 


P B (r,0,t) = (1/A B A B ) {F B (t) 6(0-6°) 6(r-rJ) 


B. 


B, 


( 65a) 


( 65b) 


+ F B ( t ) 6(0-0“) 6 ( r-r B ) } 

where T,B denote the external and internal loads, 6 is the Dirac 

B T 

delta function and for a circular plate [26] A^' = 1 , 

B T 

A 2 f = r . The generalized random forces corresponding to point 
loads given in Eos. 65a and 65b are 


,%,<»> = f?>> x sq ( r 1 ' e 1 ) 


(66a) 


B 


P sg (») = (F®( U ) x* g (rf,e*) + ?!j«o) x^ g (r®,e®)l 


(66b) 


Following the procedures of Ref. 46 and assuming the point 
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loads are stationary and independent, the spectral densities of 
normal plate deflections w T , w B can be determined from 


OO OO OO 00 


s3' B (r, Q , 

w 


- l l J • J- ^T,B 0 sq ( T,B 0 jk ) T,B S sqjk (w) (67) 
s=0 q=l j=0 k=l 

+ P A sq^P A jk^B,T S sqjk^ X sq* X jk^^sq® jk 


T , B°sq ( a) ) = { H gq B ( a) )/ m T/B }/ <5 s q ( “1 


p A sq(o) = T 0 sq (o)) { Hgg( u )/ m B } (kg + (l/6)m s w ) 


( 68 ) 


(69) 


p * sq (u) = 1 - {H sq (co)H B q U)(k s + (l/6)m sW 2 ) 2 }/m T m B 


(70) 


The asterisks in Eq. 67 denote complex conjugates and B,T s sq jj < 
are the cross spectral densities of the generalized random 
forces. 

For the two stationary independent point loads acting on the 
external plate it may be shown that the cross-spectral density of 

Q 

the generalized forces T sqjk(co) may be determined from 


T S s<ijk ( “ > = (4 < u) x sq (r l' 9 ?' ) x jk (r l' e I> + 


(71) 


S F 2 <“’ X s, (r I' 0 2 » X Jk' r 2' 0 2>) 


T T 

where S„ and S_ are the spectral densities of the point loads 

F i t 2 

T T 

Fi and F 2 . Similar expressions can be developed for point 
loads acting on the interior plate. 
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The vibrations of the bottom plate generates noise inside 
the cylindrical enclosure. By setting w B = w L at x = 0 and w B 
= w R at x = L, substituting these results into Eqs. 6 and 7 , the 
solution for interior noise pressure P 2 can be obtained. 

3.8 Response of Stiffened and Interconnected Structures 

The interconnected structural system shown in Fig. 7 is 
composed of a discretely stiffened panel and a cantilever 
stiffened box beam. The beam is attached to the stiffeners of 
the skin-stringer panel. The displacement response of the 
stiffened panel is needed for the solution of the perturbation 
pressure p given in Eq. 20. To develop a solution procedure, 
the structural system shown in Fig. 7 is separated into a 
stiffened panel, stiffened beam and a coupled panel-beam 
structure. The solution procedure is based on transfer matrix 
techniques. 


3.8.1 Response of a Stiffened Panel to Random Loads 

The vibrations of various panel units inside the 
habitability modules might be generated by various mechanical 
and/or electrical equipment that are attached to these panels. 
Consider a flat stiffened panel located on a sidewall at 
z = 0 , a o < y < a o + L , bo < x < b o + L is simply supported at 

y X 

the edges normal to the stiffeners as shown in Figs. 7 and 8. 
Assume the random loads from the substructure are transferred to 
the panel at the stiffeners through shear force, bending moment 
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and twisting moment action. Now consider a panel exposed to a 
random pressure p r (x,y,t), random point load Fg(t) and random 
couples My(t) and M x (t). The governing equation of motion for 

this single bay panel is 
DV 4 w + cw + m s w = p r (x,y,t ) + F 0 (t) 6(x — x a ) 6(y - y „ ) + 

M y (t)S(x - x 0 ) S'(y- y Q ) + M x (t) 6'(x - x 0 )S(y - y 0 ) (72 ) 

where D is the plate stiffness, c Q is the visous damping 
coefficient, m s is the panel mass per unit area, w(x,y,t) is the 
normal deflection, 6 is the Dirac delta function and the prime 
denotes a derivative. The point couples acting on a plane 
parallel to the axis are indicated by the subscript associated 
with them. Note that in formulating Eq. 72, the effect of 
acoustic radiation pressure is neglected. 

The solution for the panel deflection can be written as 

09 

w ( x , y , t ) = )’ q ( y , t ) X ( x ) (73) 

n=l n n 

where q n are the generalized coordinates and X n are the modes 
corresponding to the x-direction. Assuming simply supported 
edges at x = 0, L x , the modes are 

X (x) = sin (£2i) (74) 

x 

Taking the Fourier transform of Eqs. 72 and 73 and following 
the Galerkin procedure, the differential equation for the 
generalized coordinates q (y,aj) can be obtained [34,35,43]. 
Introducing relationships between various derivatives of 
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qn(y,w) in terms of slope, bending moment and shear, the 
solution can be written in convenient state vector form as 


{Z n } « 



( 75 ) 


where 6 , 9 , M , and V are the modal components of deflection, 
n n n n * 

slope, moment and shear, respectively. Then, the response state 
vector at station j on the panel is 


(Z 


n 




[F] . (Z } r . . 
3 n j-1 



[FCYj 


5)1 (K n (?)}d? (76) 


where the superscripts 1 and r indicate either the left or right 
of station j, respectively, [F] is the field transfer matrix 
which transfers the stste vector across the panel [33,37,46] and 

is the matrix of generalized random forces. 

Now that we can transfer the state vector across the panel, 
we next develop procedures to transfer the state vector across 
the stiffener. The stiffener does not interfere with the 
continuity of the deflections and slopes in the skin on either 
side of the line of attachment between the stiffner and the 
skin. However, the stiffner, because of its elastic and inertial 
properties, produces an abrupt change in the moment and shear of 
the skin at the line of attachment. Since the substructures are 
attached to the stiffeners as shown in Fig. 7, axial loads will 
be induced along the coordinate x. The details of developing a 
transfer matrix for the stiffener can be found in Refs. 37 and 
43. Then, the state vector at any arbitrary location s, where 
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( 77 ) 


s = y + 1 y , can be expressed as 
s m=l m 

where matrix ^ [T] o transfers the state vector from station 0 to 
s 

station s such that 

l s {T]o — [F a ] \Ga\} [F]j [Ga].?'-i * ‘ ' mi [^a]o • ( 78 ) 


in which [G A ] is the point transfer matrix for the jth stiffener, 
[F g ] is a field transfer matrix which transfers the state vector 
over a portion of a panel located between stations j and j + 1. 
Transfer matrix ^{E}o represents the effect of the distributed 
and concentrated loads acting on the panel and the axial loads 
acting on the stringers: 


[{E } 1 0 = i [T \ l 0 {L q } + l JT\\ {Li} + • • • + i[r]J {Lj} + {L s } 


where the loading matrix Lj 
{L.} = {L 9 .} + {L ^ } 


(79) 


(80) 


in which 


{L f .} = / 0 J [F(y. - ?)] (K n (C)}d5 


(81) 


P 

and matrix {L^ } includes the effect of the axial loads induced 
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in the stiffeners [43] . 

The solution for the state vector {Z }n in Eq. 77 can be 

n 

obtained by utilizing the boundary conditions at the extreme ends 
of the panel, i.e. at y = 0 and y = L, and by extending Eq. 77 to 
the length of the panel. In this approach, simple, fixed, free 
or elastic supports at the boundaries can be included. For the 
case when either simple, fixed or free supports are used at the 
ends (not necessarily the same at both ends), the response state 
vector at location s can be expressed as 


i 


= - 


3 


tlk 

hk 

tzk 

- t 4 k 


til 

t2l 

tzi 

t 4 l 




(82) 


where the t^j element is the (i,j) element of transfer matrix 
[T] ; E^ is the i-tn element of matrix E ; and (k,l) denote the 
matrix elements corresponding to boundary conditions at the 
extreme left end (station 0): 


(k,l) 


(2.4) ) 

(3.4) > 

( 1 , 2 )) 


simple 

fixed 

free 


(83) 


Similarly, (e,f) denote matrix elements according to boundary 
conditions at the extreme right end (station N): 
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((1,3)/ simple 
(e,f) =-<(1,2)} fixed 
((3,4)) free 


Finally, the response- vector in terms of displacement, 
slope, moment and shear can be obtained from 

CO 

{&(*>")}!= LWi-M*) (85) 

n = l 

The displacements from Eq. 85 can be used in Eq. 20 to obtain 
sound pressure inside a rectangular enclosure. 


3.8.2 Response of a Stiffened Beam 

Consider a discretely stiffened beam composed of piecewise 
continuous segments as shown in Fig. 9. Concentrated masses, 
translational, torsional, longitudinal springs and point dampers 
can be added to the beam structure at any arbitrary location. 
Furthermore, random loads (distributed or concentrated) can act 
at any arbitrary location on the beam. 

The deformations, forces, bending and twisting moments of 
the elastic beam at a given location s can can be represented by 
a state vector {Z^}^ • Then, the response state vector can be 
obtained as with the skin stringer panel. Thus, we write 

{Z b } l s = l s[T b ] l 0 {Z b } l 0 + l 3 {E b } l 0 ( 86 ) 

where the transfer matrices now correspond to those of a 
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where (4^/ 6^, 0^} are angles of twist about the elastic axis, 

vertical displacement, slope in the vertical plane, respectively, 

and station K is at the base of the connection between the 

cantilever beam and stiffened panel. The sub matrices 
l u ll 

{T^ } and {T^ } are (3x3) reduced transfer matrices of a beam 
with bending and twisting modes and free boundary condition at 
z=0 for which the components of the state vector are 
(Mfo, V b , } 0 = {0, 0, 0} . The second part of Eq. 87 reflects 
the effect of point loads which can be loacted at any arbitrary 
number of stations r = 0 , 1 , 2 , . . . , K-l . Similarly, the submatrices 



the input load vector {M 0 ,F 0 ,M T }r in which M 0 ,F 0 and M T are the 
point couple, point force and twisting moment acting at station 

r. 

In the formulation of the skin-stringer problem, it was 
assumed that in-plane deformations with respect to directions x 
and y are small and can be neglected. To satisfy the 
compatibility conditions between the skin-stringer and cantilever 
beam at the point of attachment K, we assume that for the beam 
1 K = 0 • Utilizing these conditions and Eq. 87, we 

obtain 


06 

<56 

06 






+ E 1', 

T = 0 



( 88 ) 


At point K where the beam joins the panel, compatibility and 
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equilibrium conditions need to be satisfied. Assume the 
cantilever beam is connected to a stiffened panel at a distance 
x = a as indicated in Fig. 7. The slope of the panel-stringer 
response along the x coordinate can be obtained by- 
differentiating the displacement component of Eq. 77 with respect 
to x. Then, 


OO 

{4>b} K — ^ K^ 11 * 12 *13 

71=1 



mxa 
cos — — 
■Lj x 


(89) 


in which t^j are elements of transfer matrix [T] corresponding to 
the skin-stringer panel, Eq. 78. Inserting Eq. 89 into Eq. 88 
and using the results in Eq. 87, the response bending moment, 
shear and torsional moment at the connection of the beam to the 
skin stringer panel can be written as 



£ ( - Ps'i [zJTW*] + p? r 'i) 



+ PVIPVT 1 


£[r' I Wig) 


nixa 


cos 


where transfer matrix [ T ^ ] is 


(90) 


[T t \ = 

K 


0 

0 

t li 


0 0 

0 0 

*12 *13 



14 


J 0 


(91) 


Consider now that the cantilever beam is connected to two 
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adjacent stringers on the panel and that all the loads from the 
beam are transferred into the stringers at stations j-1 and j. 
These loads, given in Eq. 90, are functions of {Z n >o and are 
inputs to the skin-stringer panel system. Then, the response 
state vector at station x of the stiffened panel is 

{Zn}‘,= '.Pi!, {ZnYo - M-l {Lin} ~ K { W ( 92 ) 


where 

= {0,0,0, Vj f/n } 


(93) 


in which 



nma 


s 

2 


/\~ 


nrca r 


nna 


i / "• i 


(94) 


In obtaining Eqs. 92 and 94 it was assumed that each stiffener to 
which the beam is attached shares one-half of the shear. The 
torsional moment produces loads in the z direction which can be 
approximated by four forces located at the four corners of the 
cross-section, each of which is equal to XM T [56,57], Also, 
otj, denote the distances in the x direction of the upper and 

lower corners of the cross section. 

Now for the sake of numerical tractability, we proceed to 
examine the case where response of the panel is dominated by the 
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bending moment {M^}^ . Then, Eq. 94 reduces to 

fnn\ nwa 

T “ = v ' "« = “ \m) cos 17 


which implies that 


{L n } = {Li t un} = {0,0,0,v n } 


Using Eqs. 90, 93 and 95 we can rewrite Eq. 92 as 

{£„}' = , M{z»Yo + , .\d,»] ( eWi { nl 

\ T - 0 K ,T > 

/ oo 

+ i( Dnin] I ^ [ Dr/m } {ZmYo 


m=l 


where 


iP/„] = -(;[rj;. 1 + '[r]j) {L„) 
iDnr] = -it;- j |r k 1 “)-‘[rj2»j + [J-2.J 


l s li 

in which [T, ] is a (lx 3) matrix obtained from [T, 

D D 

2 g 

eliminating the second and third rows; [T^ ] is a (3 

2 r 

obtained from [T^ 0 ] by eliminating the third column; 

r 2 1 

(1x2) matrix obtained from [T^ ] by eliminating the 

r 

third rows and third column. Also, 


(95) 

(96) 


(97) 


(98) 

(99) 
by 

2) matrix 
2 s 

[T b s ] is a 
r 

second and 
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( 100 ) 




[D]Vm 



nna 
cos — — 

L x 


[T'] 


( 101 ) 


where [T^] is given in Eq. 91. Note that the dimensions of 

matrices ^[T]o, ^ [D x ] [D TT ], [D TTT and [D TT7 ] are 
s u s In s IIt Illn IVm 

respectively, 4x4, 4x1, 1x2, 4x3 and 3x4. 

The state vector (Z >o are (z , , z .} , where k,l are given 

n nK n x 

in Eq. 83 according to the boundary conditions at the extreme 
left end of the panel. Now we approximate the infinite summation 

oo V 

£ with \ where v is a selected integer number. Then, 
n,m=l m,m=l 

using Eq. 97 for the entire panel and the boundary conditions at 
the extreme right end, a system of equations can be written in 
the following form: 


' Z\k ' 

l 

Nli 

Z 11 


N2i 

z 2k 


N 1 2 

Z 21 

> = i 

N2 2 

Z v k 


Nl v 

- Z V 1 J 

0 

v N2 u j 


( 102 ) 


where 



( 103 ) 
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( 104 ) 


/£ 1 1 -j-Hln 

Hii+5in 

R 1 1 2 

SI 12 

* * * Rx 1 1 / 

S' lii/ 

K2 l +R2 li 

H2 1 +S2 ll 

R 2 12 

S2\2 

R2 lu 

S 2 1 ,, 

HI 21 

5121 

KI 2 ~\~R I 22 

H12+S122 

• • • R%2u 

5 121 / 

.R221 

S2ai 

K 22 -\-R2i22 

H 2g -b<S 222 

. • • R22v 

S 221 / 

Rl»i 

Sl„i 

RXv2 

Slu2 

♦ . . Klv+Rlw 

H 1^ +S \ vv 

H2,/i 

S2 v i 

RIi/2 

S2 U 2 

• •• K2 tJ -{- R2 VV 

H2,/ + S2 l/l/ 


and 


•Rlnm 

R2 n m 

Slnm 

Slum 

~ N\^IIIn\ 

[D 

'Kl n 

Hl n 

^ T [ tek 

tel 

K1 n 

H2 n 

N Jfk 

tfl 


(105) 


(106) 


in which the superscripts e,f in matrices N^IIIn^ 

denote rows and k,l in matrix [Dj Vm ] indicate columns retained 
from the total matrix according to the boundary conditions at 

cf af 1 nTIQ Tvf and O rponpof 1 tto! u r PV"» T 7 a 1 1 1 qc P arud \r 1 v 

v* vlt'liO i.'i uiiu V-/ f i-CO|JCuuxvvJ.j • j. nt v oa -i- C/ O w j_ / j_ unu JS. f J. d i_ O 

given in Eqs. 84 and 83, respectively. For example, if the panel 

is simply supported at station N, the first and the third rows of 
matrix J;[D T ] are retained. Elements t..- • in Eq. 106 are defined 
as in Eq. 82. 

The natural frequencies of the coupled skin-stringer-beam 
system can be obtained by setting the determinant of matrix [W^] 
to zero, i.e. , 


W 


= A( co) = 0 


(107) 
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Finally, from Eq. 102, the unknown elements of the state 
vector at station 0 can be obtained. Then, using Eq. 97, the 
final expression for the response state vector at any location s 
on the panel can be expressed as 


{Zn}‘ 3 - i[^4n] ^ [Diir] j 


(108) 


where 


iW “ { Z Z }' + iHM + UDm.) tlDiU { £ } 


(109) 


Consider now that the bending moment of the stiffened beam 
is acting in the plane parallel to the y axis. The slope of the 
panel at the vertical plane along the y coordinate is obtained 
from 


{<Mjc 


y : k [*21 *22 (23 *24]o {ZnYo s i n 


( 110 ) 


where Eqs. 95 and 96 are now replaced by 


y 


n 



. n na 


sm 


L x 


{M b } K 


(111) 
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which implies that 


{L n > = (0, 0, v n , 0} (112) 

All other expressions remain the same, but now, instead of Eq. 
101, we have 


[DiVm 



. nna 
sin — — 

L x 



K 


<21 


0 

0 

< 22 


0 

0 

<23 



<24 


J 0 


(113) 


As a final case we considsir axial for^os act i ny on a beam 
producing longitudinal vibrations in the beam and transverse 
vibrations of the skin-stringer panel. Now, the compatibility 
conditions of the skin-stringer panel and beam interface of the 
displacements in the z direction are satisfied by setting 


nna 


{u b } K = *12 <13 < 1 4 1 0 {Zn} o Sin — 


(114) 


n= 1 


Equation 96 remains unchanged, but Eq. 94 needs to be modified to 
include axial forces. It can be assumed that (1) the axial load 
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is transferred into the panel as four point forces acting at 
corners of the cross-section, in which case Eq. 95 is replaced by 


n 7r a 


V n = 


2 L, 


sin 


V 2 ) . rnr(a + h b /2)] 

+ sin {A}, 


(115) 


where h^ is the length of the beam in x direction; or (2) that 
the axial load is transferred as distributed forces along the 
stringer where now 


mrhb 


co. + cos {Pb)i 


(116) 


The final solution takes the same form as that of Eq. 108 but Eq. 
101 is replaced by 


[DlVm, 


. mna ... 
sin — [T*] 


(117) 
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4.0 NUMERICAL RESULTS 


The numerical results were obtained for single wall 
cylindrical shells, double wall cylindrical shells, double wall 
composite shells, stiffened shells, double wall circular plates, 
rectangular enclosures with vibrating stiffened sidewall panels 
or partitions and rectangular enclosures of which stiffened 
sidewalls are interconnected to beam-like type structures. 
Damping of the shell structures was assumed to be composed of 
material damping, viscous damping caused by radiation effects, 
and structural damping of the core material (double wall 
configuration). Using the complex elastic modulus approach, we 
have 


E = E R (1 + ig) (118) 

where E is the complex modulus, E K is the real component of E , 
and g is the loss factor. The viscous damping coefficients were 
expressed in terms of modal damping ratios C e and C 1 
corresponding to the external and internal shells. Damping in 
the core is introduced through the core stiffeness constant 
k = k (1 + i g ) where g,, is the loss factor for the core 
material. Similar damping model was used for the double wall 
circular plates. For stiffened panel and stiffened beam 
vibrations, damping was introduced by replacing elastic modulus E 
and shear modulus G with E(1 + ig) and G(1 + i g) . 

Furthermore, the translational, k T j, and rotational, kpi, spring 
constants of the cantilever beam construction v/ere replaced with 
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k T ^(l + i gi) and k R ^(l + i g^) , respectively, where g-j_ is the. 
loss factor for the external springs. The effect of viscous 
damping was neglected for those structures. 


4.1 Noise Transmission Into a Cylindrical Enclosure 

4.1.1 Single Wall Aluminum Shell 

The aluminum shell shown in Fig. 3 has the following 

dimensions: L = 300 in., R = 58 in. and h = 0.1 in. Both ends 

are closed and the interior walls are lined with a layer of 

porous acoustic material [40] . The inputs to the shell are 

exterior random point loads located at X x = x 2 = 150 in. , 

oo 

0 i = -90 and 02 = 90 . These point forces are characterized by 

truncated Gaussian white noise spectral densities 


Fir 


F® 


0.84 lb /Hz 0 < f < 1000 Hz 
0 otherwise 


(119) 


The elastic modulus, Poisson's ratio and the material densities 

were respectively, E =10 . 5xl0^ps i , v = 0.3 and p m = 0.000259 

o /•_ 4. The speed of sound, c, in the enclosure, the air 
lbf-sec^/ in ’ 

density p and flow resistivity of porous acoustic material 

lining, , of the interior surfaces are taken to be c = 13540 

_7 2 4 , 

in/sec., p = 1.147 x 10 Ib^-sec /in and R^ = 3.74 x 10 J lbj- 

sec/in.^. The sound pressure levels are computed inside the 

o 

shell at x = L/2, r = 28 in. and 0 = 45 . The structural modal 

damping coefficients were taken to be constant, 

C = Co = 0.04 . The deflection response levels RL are 
mn 

calculated from 
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RL ( x / 0 , to ) 


( 120 ) 


2 

= 10 log [S w (x,e,co) Aco/w ref ] 

where S w is the deflection response spectral density and w re £ = 
h. Numerical procedures were developed to calculate modal 
frequencies, deflection response spectral densities and interior 
sound pressure levels. Pressurized and unpressurized modules 
were considered. 

The deflection response levels calculated at x = L/2 and 

o 

9 = 45 are given in Fig. 10 for pressurized and unpressurized 
shell. The effect of pressurization is to shift some of the 
lower frequency response peaks to a higher frequency values. 
Similar results are shown in Fig. 11 where sound pressure levels 
are plotted versus frequency. Since interior sound pressure is 
dominated by the modes (structural and acoustic) above 300 Hz, 
pressurization does not have much of an effect on sound pressure 
levels. It should be noted that the magnitudes of the point 
loads chosen in this study are arbitrary values. When the actual 
input force values are known, the results presented in this 
report can be scaled accordingly. For example, if the actual 
input forces are ten times less than the ones used in this study, 
the RL and the SPL levels given in Figs. 10 and 11 will be 
reduced by ten decibels. 

4.1.2 Double Wall Aluminum Shell 

Numerical results presented herein correspond to the double 
wall sandwich shell shown in Fig. 4 and input point loads as 
given in Fig. 6. The magnitude and location of the inputs acting 
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either on the exterior shell or the interior shell are assumed to 
be the same as prescribed by Eq. 119. The thicknesses of the 
external and internal shells are h E = 0.032 in and hj = 0.1 in. 

The stiffness, material density and thickness of the soft core 

O —6 2 4 

are k g = 4.17 lbf/in J , p g = 3,4 x 10 lb^ - sec /in and h s = 2 

in., respectively. 

The natural frequencies of a double wall aluminum shell are 
given in Fig. 12. For the double wall shell construction, the 
flexural (in phase) and the dilatational (out of phase) modes are 
included. The highest modal frequency is that of "breathing" 
mode for which n = 0 and m = 1. Results plotted in Fig. 12 
indicate that for the large shell dimensions and the ratio 
radius/length = 0.1933 chosen in this study the modal frequencies 
at n = 0 seem to converge to a single point for all values of m = 
1,2,..., 10. This suggests that in the vicinity of the 
"breathing" mode frequency large numbers of structural modes 
could couple to acoustic modes resulting in high levels of noise 
transmission. The deflection response levels of the outer shell 
are given for pressurized and unpressurized conditions in Fig. 

13. These results correspond to x = L/2, 0 = 45° , g s = 0.05, 

Co = 0.04 and inputs provided by two point loads acting on the 
interior shell at x\ = x\ = 150 in., ©i = -90° and 02 = 90° . 

As can be seen from these results, shell response is dominated by 
low frequency modes in the range 0-50 Hz. In Fig. 14, sound 
pressure levels in the cylindrical enclosure generated by two 
interior point loads are given for reverberant (hard interior 
walls) and absorbent (interior walls treated with soft acoustic 


45 



materials) conditions. For reverberant conditions, noise levels 
inside the cylinder become relatively large and are dominated by 
peaks at acoustic resonant frequencies. The effect of 
pressurization on noise transmission in the interior is shown in 
Fig. 15 for absorbent interior wall conditions. Those results 
tend to indicate that pressurization plays only a minor role on 
noise generation inside the shell. It is expected that the 
habitability modules will be pressurized for which static 
pressure differential Ap = 11.8 psi . The interior sound 
pressure generated by identical inputs but acting at different 
locations is given in Fig. 16. As can be observed from these 
results, noise generated by shell vibrations is a function of 
location of random point loads. 

4.1.3 Double Wall Composite Shell 

The dimensions of the double wall composite shell are the 
same as those of double wall aluminum shell. The outer shell 
consists of three laminae while the inner shell is composed of 
ten laminae. Fiberglass and graphite fibers are used to 
reinforce the plexiglass material. The ratio of fibers volume to 
the plexiglass volume is 0.2. The fiber orientation is 
prescribed by angle a as shown in Fig. 4. The elastic moduli, 
Poisson's ratios and material densities are = 7.75 x 10^ psi, 

= 0.33, p f = 0.0002 lbf-sec^/in 4 , E g = 10.5 x 10 7 psi, 

v = 0.33, p = 0.00015 lb f -sec 2 /in 4 , E = 2.35 x 10 5 psi, 
y y ir 

v = 0.35, p = 0.00011 lb^-sec^/in^ , where the subscripts f,g,p 
P P 

denote fiberglass, graphite and plexiglass, respectively. The 
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fiber reinforcement — same pattern is used for internal and 
external shell — is arranged as follows: 1st layer fiberglass, 2nd 
layer graphite, 3rd layer fiberglass and so on. The fiber 
orientation for the three laminae of the exterior shell are 
a = -45°, 45°, -45° (first case) and a = 90°, 0°, 90° (second 
case). The fiber orientation for the ten laminae of the interior 
shell is arranged in an alternating order with 
a = -4 5 0 , 45 0 , -4 5 0 , 4 5 0 , etc. (first case) and 
a = 90°, 0°, 90°, 0° , etc., (second case). 

The natural frequencies of double wall composite shell are 
plotted in Fig. 17. A comparison of modal frequencies of 
aluminum and composite shells shows that depending on fiber 
reinforcement orientation, significantly higher modal frequencies 
can be obtained for a composite shell. However, the mass of the 
composite shell is about 50% less than that of the aluminum shell 
while all other geometric parameters remain the same. The 
deflection response levels for point loads acting on the exterior 
surface are shown in Fig. 18a for modal damping coefficient 
Co = 0.01 and loss factor in the core g„ = 0.02. As can be seen 

o 

from these results, a large number of flexural and dilatational 
modes are excited by point loads. Due to the large number of 
participating modes and modal frequency overlaps as shown in Fig. 
17, it is difficult to identify the response peaks corresponding 
to dilatational frequencies. However, for n = 0 the flexural and 
dilatational frequencies are well separated. A direct comparison 
of these results to the results given in Fig. 13 indicate that at 
most frequencies the response levels of the composite shell are 
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lower than when compared to the response levels of the aluminum 
shell. However, at some frequency values the opposite is true. 
Similar results are presented in Fig. 18b but for the point loads 
acting on the interior shell. As can be seen from Figs. 17 and 

18b, response levels at the first three peaks are about the same 

for both of these cases. However, significantly different 
vibration levels might be observed at other frequencies when the 
input point loads are moved from external to internal shells. 

The location and magnitude of these loads are the same for both 
cases. 

To demonstrate the effect of shell and core damping, results 
are presented in Figs. 19 and 20 for constant modal damping 

ratios C E n = C* n = C Q = 0.04 and g g = 0.1 . The point loads are 

acting on the interior shell for both of these cases. By 
increasing modal damping of the interior and exterior shell from 
0.01 to 0.04, about 12 dB of response reduction can be gained at 
most modal frequencies. As can be seen from Figs. 19 and 20, 
only about 2-4 dB of the response reduction is achieved at some 
peaks when damping in the core is increased from 0.02 to 0.1. 
However, the shells forming a double wall construction are bonded 
to the core. Thus, the cumulative effect of damping on vibration 
response would be similar to the combined results given in Figs. 
19 and 20 . 

Figure 21 depicts sound pressure levels for an aluminum and 

fiber reinforced laminated shell under exterior point load 

inputs, with C E = C 1 = 0.01, g = .02 and 
mn mn 3 s 

_ 8 2 

P = 1. x 10 rad-sec/in . As can be observed from these 
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results, the noise levels generated by a composite shell are 

higher than the noise levels for an aluminum shell at most 

frequencies. The mass of the composite shell is about one half 

of the mass of the aluminum shell. However, the composite shell 

is much stiffer than the aluminum one. For a shell structure, a 

shift in modal frequency could induce different coupling between 

structural and acoustic modes. The effect of structural and 

acoustic damping on sound generation is illustrated in Figure 

_ 8 2 

22. These results correspond to 3 = 1 x 10 rad-sec/in . As 
can be seen from these results, a significant amount of noise 
reduction can be achieved in a composite shell by increasing 
structural and acoustic damping. The results shown in Fig. 22 
indicate that for acoustically hard interior walls (Z a -* ®) , 
the noise levels in the cylinder become relatively large. 

A direct comparison of interior sound pressure levels in the 
cylinder excited by exterior and interior point loads is given in 
Fig. 23. The loading conditions are the same for both cases. 
Since vibration coupling is provided by the viscoelastic core, 
the noise generated in the interior is a function of how the 
point loads are acting on the double wall shell. The results 
presented in Fig. 23 correspond to point loads acting on the 
interior shell at x\ = ~ L/2, 9 i = -90° and 02 = 90° . The 

fiber orientation of the three layers (Fig. 4) at the exterior 
shell is described in Fig. 24. The fiber orientation for the ten 
layers of the interior shell are: (A) 0 ° , 22 . 5 0 , 45 0 , 45 ° , 22 . 5 ° , 

0°, 90°, 90°, 90°, 90° (B) 90 ° , 0 ° , 90 ° , 0 ° , 90 ° , 0 ° , 90 ° , 0 ° , 90 ° , 0 0 (C) 
-45°, 45°, -45°, 45°, -45°, 45°, -45°, 45°, -45°, 45°. These results show 
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that shell response and interior noise are functions of fiber 
orientation in a composite shell. The interior noise levels 
might be tailored to meet specific needs by selecting a suitable 
fiber orientation. However, interior noise is a function of 
frequency and only specific frequency bands might be affected by 
this procedure. 


4.1.4 Discretely Stiffened Cylindrical Shell 

Analytical formulations and noise transmission calculations 
have been performed for the case where the stiffened shell shown 
in Fig. 5 is represented by an equivalent orthotropic shell 


model. The solution procedure for shell response and noise 
transmission is similar to a monocoque shell presented in Secs. 
3.1, 3.3 and 4.1.1, but the natural vibration frequencies and 
mass of the stiffened shell structure are calculated by the 


methods presented in Ref. 10. 
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frequency (below 1000 Hz) vibrations and noise generation, such a 


model might be adequate to evaluate noise transmission 


characteristics of space station habitability modules. The 
structural parameters chosen in this study are typical of the 
proposed habitability modules where L = 420 in., R = 78 in and h 
= 0.1 in. Skin and stiffening elements (rings and stringers) are 
constructed from aluminum with the following section and material 
properties : 

A r = (cross-sectional area, ring) = 1.897 in^ 

2 

A g = (cross-sectional area, stringer) = 0.252 in 

I r = (moment of inertia of ring about its 
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centroid) = 5.294 in^ 

I = (moment of inertia of stiffener about its 

o 

centroid) = 0.255 in^ 

J r = (torsion constant for ring) = 0.1152 in^ 

J = (torsion constant for stiffener) = 0.000302in^ 

fa 

E = E r = E s (moduli of elasticity) = 10.0 x 10^ psi 

G = G„ = G„ (shear moduli) = 3.846 x 10^ psi 

L fa 

v = v = v = (Poisson's ratio) = 0.3 
r s 

p = p = p = (material density) = 0.000259 lb- 
s r s 

2 /• 4 

sec vm 

The natural frequencies of the orthotropic shell are shown 
in Figs. 25 and 26 for two different cases of structural 
configuration. The modal frequencies for the first few 
circumferential modes (N = 0,1,2 ,3) are not affected much by the 
number of ring frames. However, for mode numbers larger than N 
=3, the ring frames have a strong effect on modal frequencies. 
Furthermore, for the circumferential wave numbers larger than six 
and longitudinal modes higher than ten, the modal frequencies 
tend to converge to a single line. From these results, it can be 
seen that at each selected frequency several modes could be 
contributing to structural response and noise transmission. The 
noise transmission inside the cylindrical enclosure shown in Fig. 
5 was calculated for a variety of structural configurations. The 
inputs are point loads acting at xj = X 2 = L/2, 

0! = 90°, 0 2 = -90° . These inputs are assumed to be stationary 
Gaussian white noise random processes characterized by truncated 
spectral densities 
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0.01 lb 2 /Hz 0 < f « 1000 Hz 
0 otherwise 


( 121 ) 


It is assumed that some absorbent acoustic material is intact at 
the interior walls of the cylindrical enclosure. 

The generated interior SPL in the enclosure at x = L/2, 

9 = 45° and r = 63 in are presented in Fig. 27 for several 
geometric stiffening configurations. For the stiffened cases it 
is assumed that only transverse rings are present as stiffening 
elements. As can be seen from these results, at low frequencies 
(below 150 Hz) interior noise levels for unstiffened shells are 
higher when compared to the results of a stiffened case. 

However, for frequencies above 150 Hz, higher noise levels might 
be generated at some frequencies for a shell stiffened with 
rings. This could be due to the fact that the acoustic modes at 
these frequencies are strongly coupled to the shell structural 
modes. The sound pressure levels for a pressurized shell 
stiffened with 10 rings and stringers which are spaced at 10 
inches apart are given in Figure 28. The results tend to 
indicate that more interior noise is generated by a thinner 
shell. The effect of stringers on interior noise levels is 
illustrated in Fig. 29. As can be seen from these results, the 
sound pressure levels do not change by much when stringers are 
added to the positions shown in Fig. 30. These results 
correspond to a shell stiffened with ten heavy frames and 
stringers which are spaced 10 inches apart. The sensitivity of 
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noise generation due to different locations of input forces is 
clearly demonstrated in this figure. Structural modes which are 
not excited for a particular forcing condition might become 
efficient sound radiators for a different set-up of input forces. 

Numerical results were also obtained for a shell stiffened 
with relatively small ring frames. The structural parameters 
selected are: 

A r = 0.228 in 2 
A s = 0.295 in 2 
I r = 0.4437 in 4 
I s = 0.5087 in 4 
J r = 0.000122 in 4 
J s = 0.000246 in 4 

The material properties are the same as given in the previous 
examples. A comparison of interior sound pressure levels for the 
two cases of different ring frame stiffening is shown in Fig. 

31. These results indicate that interior noise levels are 
significantly higher for a shell stiffened with small frames 
(frequency range 50-500 Hz). The results shown in Fig. 32 show 
the effect of increasing the number of small frames from 10 to 
41. In the frequency ranges 50-200 Hz, 400-500 Hz, interior 
noise levels are higher for a shell stiffened with 41 small ring 
frames than for a shell stiffened with 10 heavy frames. However, 
for frequencies above 750 Hz more noise is transmitted when ring 
frames are large. 
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The results presented indicate the sensitivity of interior 
noise environment inside the habitability modules due to changes 
in structural, geometric and loading conditions. These results 
were obtained for a particular level of the point load 
intensity. At the present time, the magnitude characteristics 
and location of mechanical inputs that will be present in the 
habitability modules during orbital operations are not known. 

The measurements obtained for the Skylab operations shown in 
Figs. 33 and 34 indicate typical interior noise levels. In 
general noise levels generated by various mechanical or 
electrical components are relatively low, but the total level 
might reach about 80 dB at some frequencies. 

4.1.5 End Plates 

Consider noise is being generated by vibration of double 
wall end plates which are excited by point loads as illustrated 
in Fig. 6. The core separating the double wall end plate 
construction is taken to be relatively soft in order to allow for 
dilatational modes to be present. The coupled modal frequencies 
of the double wall aluminum caps for s = 0,1, 2, 3 (number of nodal 
diameters) and q = 1,2,..., 10 (number of nodal circles) were 
calculated. The first three structural modes for zero number of 
diametrical nodes are shown in Fig. 35. 

In Fig. 36, the first three radial acoustic modes are 
illustrated. These results were obtained from Eq. 14 for zero 
number of nodal diameters (i.e., j = 0, no variation in 0 - 
direction) and k = 1,2,3, where k represents the number of nodal 
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circles (r-direction ) . 

The sound pressure levels at x = L/2, r = 23 in. and 9 = 

45° due to noise transmitted through the double wall circular end 
plates located at x = L, are shown in Fig. 37 for reverberant and 
absorbent interiors. The input is a uniform 120 dB acoustic 
pressure acting on the exterior end plate. In this case, the end 
plate located at x = 0 is assumed to be rigid. The reverberant 
and absorbent conditions are simulated by selecting Z^, 

_ 7 

Z R + ®, p = 0 , and Z L , Z R as given in Eq. 9 and p = 1 x 10 

rad-sec/in , respectively. As can be seen from Fig. 37, a large 

number of acoustic modes are excited by the vibration of the end 

plates for reverberant conditions. Modal plate damping is taken 

T B 

to be constant and equal to C = C = 0.06 . The structural 

^ sq sq 

loss factor of the core g^ = 0.02. The noise transmission of 

s 

the end caps is predominantly low frequency. The fundamental 
circular plate frequency is 3.73 Hz while the lowest acoustic 
modal frequency in the shell enclosure is 22.56 Hz. From these 
results it can be seen that under uniform random pressure input, 
the noise transmitted by the double wall shell and circular end 
plates could be relatively large over the selected frequency 
range . 

The results presented in Fig. 38 illustrate the difference 
between the noise transmitted due to a uniformly distributed 
acoustic pressure input and sound generated by point loads. In 
both cases, the random excitations are acting on the exterior 
plate of the double wall construction located at x = L and the 
end plate located at x = 0 is assumed to be rigid. The uniform 
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input is 120 dB acoustic pressure and the two point loads are 
characterized by a truncated Gaussian white noise spectral 
density given by Eq. 119 and are located at rj = r£ = 28 in. and 

m rp 

01 = -90°, 02 = 90° . The absorbent conditions are described in 

Eq. 9 and the equivalent acoustic damping parameter is 
_7 9 

P = 1 x 10 rad-sec/in^ . The modal damping ratios are taken to 
be constant and equal to 0.06. The loss factor of the core 
g^ = 0.02 . The sound pressure levels are calculated at x = 

L/2, r = 23 in., 9 = 45° . From Fig. 38 it can be seen that the 
uniform acoustic pressure tends to generate more noise in the low 
frequency region while the sound generated by point loads inside 
the enclosure is about 10-15 dB higher in the high frequency 
region. 


4.1.6 Total Interior Noise 

Due to the assumption of independently vibrating double wall 
shell and end plate systems, the total interior pressure can be 
calculated by a superposition of the individual contributions. 

In Fig. 39 results are shown of noise generated inside the 
enclosure due to uniform random pressure applied on the exterior 
surfaces of the double wall shell and double wall end plates. It 
can be seen that transmitted noise is dominated by end plate 
vibrations for frequencies up to 200 Hz and by shell vibrations 
for frequencies above 200 Hz. Then, the total interior pressure 
is presented in Fig. 40. These results indicate that neglecting 
noise transmitted by the end caps would underestimate interior 
sound pressure levels for the low frequency region. Similar 
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results are presented in Figs. 41 and 42 but for random point 
load inputs. These loads were applied on the exterior surfaces 
of the shell and end plate systems. As can be observed from 
these results, low frequency noise is dominated by end plate 
motions . 


4.2 Structural Response and Noise Transmission of Rectangular 
Enclosures 

Numerical results have been obtained for several simplified 
versions of the interconnected structure shown in Fig. 7. 

Examples of such simplified structural systems are given in Figs. 
43,44 and 45. The stiffened panel is taken to be composed of 
three equal bays. The supports at y = 0 and y = are assumed 
to be elastic stiffeners. The stiffeners are identical and 
placed at equal distances, y^ = y 2 = ¥3 = 8.2 in. and L x = 20 
in. The physical parameters for the aluminum panel and 
stiffeners are typical of a transport jet aircraft [ 9 ] . For 
the special case shown in Fig. 43, the box beam is not attached 
to the skin-stringer system and inputs to the structure are 
provided by four random point forces or four point couples. The 
geometric parameters of the box beam shown in Fig. 44 are = 60 
in., A = 12.48 in^ (cross-sectional area), Iy. = 35.127 in^_ 

(second moment of area about y-axis and I v = 35.127 in^ (second 

a 

moment of area about x-axis). To reflect the effect of different 
box beam geometries, cross-sectional areas and second moments of 
inertia were adjusted for different values of beam height H and 
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beam thickness t-^ and t 2 . The sectional properties for the 
piecewise continuous box beam shown in Fig. 45 are A 3 = 5.814 
in 2 , A 2 = 9.196 in 2 , A 3 = 12.48 in 2 , I yl = 2.115 in 4 , 

I y 2 = 12.265 in 4 , I y3 = 35.127 in 4 and the distances between the 
continuous segments are z 3 = z 2 = z 3 = 20 in. Numerical results 
were also obtained for cases with concentrated mass (m 3 = 0.5 
lbf-sec /in ) , and a translational spring (k T3 = 9000 lb/in) and 

n 

rotational spring (k R i - 2 x 10 'Ib-in/rad) attached at station z 
= z 3 . Inputs to the beams are the concentrated couples M ox or 
M Qy which act at the centroid of the cross-sectional area at z = 
0. All the structural elements are assumed to be made from the 
same aluminum material as the skin-stringer panel. Damping in 
the skin-stringer-beam structure is introduced by replacing the 
elastic modulus E and shear modulus G with E(1 + ig ) and G(1 + 
ig) where g is a loss factor. Furthermore, the translational, 
krp 3 , and rotational, k R3 , spring constants are replaced with 
k>ri ( 1 + igi ) and k Rl (l + igi) , respectively where g^ is the 
loss factor for the two external springs. The effect of viscous 
damping is neglected by setting the viscous damping coefficient 
to zero. 

Numerical results were also obtained for a large rectangular 
enclosure of which the dimensions could be representative of 
interior dimensions of the space station construction. The 
geometry of such an enclosure is shown in Fig. 46. The noise in 
the interior is generated by vibrations of panels or partitions 
which can be located at any arbitrary position on the 
enclosure. Vibrations to these panels might be induced by 
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various mechanical and electrical equipment such as fans, life 
support systems, printers, experimental devices, etc. 


4.2.1 Response of Stiffened Panel to Point Forces and Point 
Couples 

The four point forces or four point couples are assumed to 
be of equal strength and characterized by truncated Gaussian 
white noise spectral densities 

(o . 84 lb 2 /Hz 0 < f < f / 

S p =} U C (122) 

i / 0 otherwise \ 

0 . 84 ( lb-in) 2 /Hz 0 < f < f 

' u 

0 otherwise 




where f is the frequency in Hz, i = 1,2, 3, 4 and f u is the upper 
cut-off frequency. These point loads are located at the 
positions as indicated in Fig. 43. The response spectral 
densities corresponding to point force inputs are given in Fig. 
47. These results are for an upper cut-off frequency f u = 800 
Hz. The panel deflection response is calculated at x = 10 in. , y 
= 4.1 in. and y = 8.2 in. The first point is at the middle of 
the first panel bay while the second point is at the middle of 
the first stringer. Depending on the value of frequency, 
significant differences in response values are obtained at 
different locations of the panel. Similar results are presented 
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in Fig. 48 for x = 5 in., y = 12.1 in. and several values of the 
structural damping coefficient g. A large amount of response 
reduction can be achieved by increasing the structural damping. 

Deflection response spectral densities for point couple 
inputs are shown in Fig. 49 for two different values of input 
conditions. From the results presented in Figs. 47 and 49, 
differences can be seen in the skin-stringer panel response 
between the point force and the point couple inputs. 

4.2.2 Response of Stiffened Panel to Calculated Bending Moment 
Inputs 

The response of a stiffened panel to calculated point couple 
inputs was obtained. These inputs were determined from the 
response solution of the cantilever beam shown in Fig. 44. The 
beam is assumed to be clamped at z = and free at z = 0. The 
bending moment spectral densities were calculated at the root of 
the cantilever beam. The spectral densities shown in Fig. 50 for 
two different cross-sectional areas of the box beam are then used 
as point couple inputs to the stiffened panel. The response 
spectral densities of the skin-stringer panel calculated at x = 5 
in. and y = 12.3 in. are shown in Fig. 51 for two different 
geometries of the cantilever beam. The inputs to the panel are 
assumed to be acting on two adjacent stringers at a distance 
a = 6.666 in. as shown in Fig. 7. The solid line given in 
Fig. 51 corresponds to a beam with a fundamental bending 
frequency of 52 Hz, while the conditions represented by a dashed 
line are for a beam with a natural frequency of 123 Hz. Since 
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the natural frequency of one of the skin-stringer panel modes is 
approximately 123 Hz, a very strong peak appears in the response 
spectral density when the beam resonant frequency coincides with 
one of the skin-stringer panel's resonant frequencies. 

Therefore, the dynamic conditions of such substructures could 
play a significant role in controlling the response levels of a 
stiffened sidewall. 

The effect of direction of the input point couple action is 
illustrated in Fig. 52. The dashed line illustrates skin- 
stringer panel response at point x = 5 in. and y = 12.3 in for 
Mqx input (bending effect on the stringers) and the solid line is 
for M Q y (twisting effect on the stringers) input. These results 
indicate that the deflection response is significantly smaller 
for bending moment input M QX than for twisting moment input 
M Q y. This is due to the fact theat these stiffeners provide more 
resistance in bending than in torsion.- 

4.2.3 Response of a Coupled Skin-Stringer-Beam Structure 

Skin-stringer panel response was calculated for the 
geometries shown in Figs. 44 and 45. The inputs are point 
couples M QX , for which the spectral densities are given in Eq. 
123, with an upper cut-off frequency of 320 Hz. A direct 
comparison of the results is given in Fig. 53 for the cases of a 
completely coupled problem (Fig. 44) and under special conditions 
where the inputs to the skin-stringer panel are the bending 
moments of a clamped-free beam. When the box beam is attached to 
an elastic stringer, the mechanism for transmitting vibrational 
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energy into the stiffened panel is different for separately 
calculating the bending moment of the cantilever beam and then 
using the results as inputs to the skin-stringer panel. Results 
indicate the response levels of the coupled problem are lower at 
most of the frequencies. 

The effect of different beam geometries on skin-stringer 
panel response is illustrated in Fig. 54 for the structure shown 
in Fig. 44. The deflection response spectral densities of the 
skin-stringer-beam structure shown in Fig. 45 are presented in 
Fig. 55 for damping coefficients g = 0.02 and g 1 = 0.05. These 
results illustrate the effect of the concentrated mass and 
elastic springs attached to the beam at station z = Skin- 

stringer panel vibrations are particularly sensitive to changes 
in beam dynamic characteristics in the frequency range of 
approximately 10 Hz to 80 Hz. The first beam bending mode occurs 
in this frequency range and moves to lower frequencies with the 
addition of mass and to higher frequencies with the attachment of 
elastic springs. By adjusting the structural and geometric 
parameters of the box beam structure, the response 
characteristics of the skin-stringer panel can be tailored to 
prescribed conditions. 

4.2.4 Noise Transmission 

For the calculation of noise transmission into an enclosure, 
the simplified models shown in Figs. 56 and 57 have been 
chosen. A cantilever box beam is attached to two stiffeners of a 
discretely stiffened sidewall. A random poit couple acts at the 
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free end of the cantilever beam. The walls at z = 0 and z = d of 
the acoustic enclosure are treated with uniformly applied 
absorptive materials which are represented by a point impedance 
model as given in Eq. 9. The results obtained here are for = 

4 x 10^ mks rayals/m. The point couple M_ v (t) acting on the 
cantilever beam is assumed to be characterized by a truncated 
Gaussian white noise spectral density given in Eq. 123. The 
upper cut-off frequency was chosen to be f u = 600 Hz. 

Numerical results were obtained for a = 142 in., b = 50 in., 

d = 48 in., a Q = 60 in., b Q = 15 in., y^ = y 2 = Y 3 = 

8.2 in., L x = 20 in., Ly. = 24.6 in., Iy = 38.127 in^ (second 

moment of area of box beam), A = 12.48 in 2 (cross-sectional area 

of the box beam), L = 60 in. (beam length) and a = 6.666 in. 

Interior noise in the enclosure was calculated at x = 25 in., y = 

71 in., and z = 24 in. Noise transmission calculations were 

obtained for a loss factor for the panel g = 0.02. The acoustic 

damping coefficient p was related to the acoustic modal damping 

ratios £ . . by 
13 * 

2 

r = J3C 

ij 2 “ijo (124) 


where Pc 2 - 2w q 5 q , in which aj Q is the lowest acoustic modal 
frequency chosen from wool wqio an <3 (i.e., co 0 1 0 i- n our 

If 

case). The damping ratio corresponding to the fundamental 
acustic mode, = 0.03 , takes into account the contributions 
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of all the damping effects of the acoustic space. Air density 
and speed of sound in the enclosure are p = 1.147 x 10 ”^ lbf- 
sec 2 /in 2 and c = 13540 in/sec. , respectively. All the 
calculations are based on a frequency bandwidth A f = 1 Hz. 

Sound pressure levels are given in Fig. 58 for an acoustic 
enclosure with a small amount of absorption at the walls 
(Z ^ •* <» ), and for the case where the walls are treated with 
porous materials for which the point impedance is given in Eq. 

9. In this case, as shown in Fig. 58, the inputs to the 
stiffened panel are four point forces acting on the two 
intermediate stringers for which the spectral densities are 
assumed to be the same and equal to 0.84 x 10“ 2 lb 2 /Hz. As can be 
observed from these results, the peaks of the acoustic modes are 
suppressed when the acoustic absorption in the interior is large. 

The results presented in Fig. 59 are for the beam-skin- 
stringer geometry shown in Fig. 56 and acoustic absorption as 
prescribed by Eq. 9. The two cases correspond to two different 
dynamic characteristics of the beam where f i is the fundamental 
bending frequency of a clamped-free cantilever beam. In this 
case, the bending moment response spectral density of a clamped- 
free beam was calculated first. Then, the bending moment 
spectral density at the root was used as the input to the skin- 
stringer panel. These results indicate that when the natural 
frequency of the beam coincides with one of the natural 
frequencies of the stiffened panel, a large amount of noise can 
be generated inside the enclosure. To minimize this, the modal 
frequencies of the substructure should not have values which are 
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close to modal frequencies of the noise transmitting sidewall. 

A direct comparison of the results of the cases of a coupled 
problem and under special conditions where the inputs to the 
stiffened panel are the bending moments of a clamped-free beam is 
given in Fig. 60. The sound pressure levels of the coupled 
problem are lower at most frequencies. These results clearly 
illustrate the significance of a complete solution when the 
elastic coupling between the substructure and stiffened panel is 
included. 

4.2.5 Noise Transmission - A Parametric Study 

Parametric studies of noise transmission were performed for 
the acoustic enclosure shown in Fig. 46. The inputs to the 
panels or partitions are random point forces characterized by an 
idealized spectral density shown in Fig. 61. This spectral 
density is composed of low level wide band noise and periodic 
tones with a fundamental tone at 60 Hz. These types of 
excitations are often produced from operations of electrical and 
mechanical devices such as fans, compressors, turbines, rotating 
and reciprocating tools, etc. Structural damping of the panels, 
acoustic damping and acoustic absorption in the interior are 
taken to be the same as in Sec. 4.2.4. Interior noise levels are 
calculated at x = 42 in., y = 210 in. and z = 30 in. (middle of 
the enclosure). 

The narrow band sound pressure levels are given in Fig. 

62. The overall and the overall A-weighted levels are also 
indicated in this figure. Similar results are shown in Fig. 63 
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where the sound pressure levels are plotted on the one-third 
octave scale. The results given in Figs. 62 and 63 correspond to 
vibrating stiffened panels located at z = 0, x = 32 in. and y = 
198 in. The panel dimensions are Yi = Y2 ~ Y 3 = 8,2 i n */ L x = 20 
in. and = 24.6 in. It is assumed that noise is generated only 
by this panel. As can be observed from these results, strong 
peaks are observed at frequencies where the tone inputs coincide 
with structural and/or acoustic resonance frequencies. Interior 
noise is dominated by a peak at 500 Hz. The results shown in 
Figs. 64 and 65 are for the same conditions as in Figs. 63 and 64 
but the vibrating panel is located at z = 0, x = 10 in. and y = 

10 in. Similar results are given in Figs. 66 and 67 for z = 0, x 
= 32 in., y = 10 in., and in Figs. 68 and 69 for z = 20 in., x = 
84 in. and y = 198 in. As can be observed from these results, 
the spatial location of the vibrating source could have a 
significant effect on interior noise levels. 

The results presented in Figs. 70 and 71 are for an 
enclosure separated by a partition at y = 140 in. (Fig. 46). 

Noise is generated by a vibrating panel located at z = 42.3 in, x 
= 32 in., and y = 0 (left end of enclosure). Sound pressure 
levels are calculated at x = 30 in. , y = 42 in. and z = 20 in. A 
direct comparison of these results with the results presented in 
Figs. 62-69 indicate that many more acoustic modes are excited in 
a smaller enclosure. 

To illustrate the effect on interior noise by vibrating 
stiffened panels of different sizes, noise transmission 
calculations were obtained for a variety of panel sizes. All 
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these panels were located at z = 0, x = 32 in. and y = 198 in. 
for the enclosure shown in Fig. 46. Noise is calculated at the 
middle of the interior space. The one-third octave sound 
pressure levels for several vibrating panel sizes are given in 
Figs. 72-74. These results indicate that the size of the 
vibrating surface could have a significant effect on the levels 
of interior noise. However, there seems no strong indication 
that a larger panel would produce more noise. Interior noise is 
controlled by coincidences of tone frequencies and modal 
frequencies of the stiffened panels. 

Noise generated by a large unstiffened panel with dimensions 
L x = 20 in., Ly = 16.4 in. is given in Figs 75 and 76. Inputs 
are two point forces acting on the panel as shown in Fig. 75. 
Since the panel fundamental frequency is at about 35 Hz, large 
peak is seen at this frequency. Similar results are shown in 


Fiyt>. 77 and 78 for the same geometric and input conditions but 


the panel is stiffened with stringers at the boundaries and at 


the middle of the panel. The effects of stiffening are lower 
noise levels in the frequency region of 0 - 300 Hz. 
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5.0 CONCLUSIONS 


Analytical models were developed to predict vibration 
response and noise transmission of cylindrical and rectangular 
enclosures to random point inputs. These enclosures are intended 
to represent the exterior and interior geometries of the 
habitability modules of the space station design. The main 
emphasis of this study was on performing various parametric 
studies of structureborne noise generation and transmission. 
Numerical results are presented for several possible input 
conditions that might arise during prolonged orbital 
operations. When the actual input conditions are known, the 
results presented in this report can be scaled to these 
conditions. The intent of this work is to give preliminary 
guidelines for constructing analytical models and evaluating 
vibration and noise levels. 

Results indicate that the shell response is strongly 
dependent on damping characteristics of the shell material and 
the core, location of the point load action, and reinforcing 
fiber orientation of the different laminae. In general, the 
response levels for a composite double wall shell are lower at 
most frequencies than those of an equivalent aluminum shell. The 
vibration response of the end caps (circular plates) are 
predominantly low frequency with the largest peak occurring at 
the fundamental mode. 

The interior noise in a cylindrical enclosure is strongly 
dependent on damping characteristics of the shell and the core, 
location of the point load action, fiber orientation of the 
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different laminae and wall absorption of the interior walls. A 
fiber reinforced composite double wall shell tends to generate 
more noise than an equivalent aluminum shell. This is due to the 
fact that the mass of the composite shell is about one half of 
the mass of the aluminum shell and increase of the modal 
frequencies of the stiffer composite shell could induce different 
coupling of the structural-acoustic modes. The noise transmitted 
by the end caps is predominantly low frequency. Thus, neglecting 
noise transmitted by the end caps could underestimate interior 
sound pressure levels for the low frequency region. Furthermore, 
by a proper selection of structural damping, reinforcing fiber 
orientation, acoustic absorption and core stiffness, a 
significant amount of lower response and higher noise attenuation 
can be achieved by a design consisting of double wall laminated 
fiber reinforced composite shells and a soft viscoelastic core. 

The response and noise transmission characteristics of a 
shell stiffened with rings and stiffeners is strongly dependent 
on the type and the number of stiffening rings. Presence of 
small longitudinal stiffeners do not seem to have much effect on 
noise transmission. Furthermore, the location of point input 
forces could have a significant effect on generated noise levels 
in the cylindrical interior. 

Transfer matrix procedures were developed to study the 
dynamic response, noise generation and transmission of stiffened 
and interconnected structures to random loads. It has been 
demonstrated that the formulation can be applied to a variety of 
discretely stiffened structures. In addition, it has been shown 
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that structural response and transmitted noise levels in the 
interior are sensitive to the dynamic characteristics not only of 
the stiffened sidewall but also of the substructure to which 
random loads are applied. ■ Significant differences were found in 
the structural response and noise transmission characteristics of 
different types of matching boundary conditions at the points of 
structural interconnection between the substructure and the 
stiffened sidewall, i.e., between a completely coupled case and 
an idealized condition where the response of the substructure is 
computed independently from the response of the noise 
transmitting sidewall. The results indicate that by tailoring 
the geometric and material characteristics of structural 
subcomponents, the vibration levels and noise transmission can be 
reduced. 

Noise generated by vibration of various interior panels 
and/or partitions is sensitive to geometric conditions of the 
interior, panel sizes, location of the vibrating panels and types 
of inputs. No simple rules seem to exist in relating these 
conditions to the levels of interior noise. The coincidence of 
structural resonance frequencies with one of the input tones 
could result in high interior noise levels. For broad bond input 
types, increase in structural damping and interior acoustic 
absorption seem to be the most effective means for noise 
control. An alternative procedure would be to isolate the 
vibrating equipment so that force inputs to the panels are 
reduced. The results presented in this study demonstrate that a 
relatively small amount of vibrational energy is needed to 
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produce relatively high structureborne noise levels 
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Fig. 1 Habitability modules of space station 


76 



owami page fs 

OF POOR QUALITY 


Standard Cylindrical Section c"o Sure 

Insulation (Typical)^ (Typical) 


Gasket 



Dome Assembly 


Specialized 

Interior 


Protective 
Shroud (Typical) 



Standard Floor 

Subsystem 
Installations 


Radiator 

Assembly 

(Typical) 




Fig. 2 Proposed geometry of habitability modules 
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Fig. 3 Geometry of a cylindrical shell model 







Fig. 7 Geometry of an acoustic enclosure and a stiffened structure. 
















PRESSURIZATION EFFECT STUDY 
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Fig. 10 Defelction levels for pressurized and 
unpressurized aluminum shells 






Fig. 12 Natural frequencies for double wall aluminum shell 







Fig. 14 Sound pressure levels of a double wall 
aluminum shell (exterior point loads) 




Fig. 15 Sound pressure levels in a pressurized 

and unpressurized double wall aluminum shell 



INPUT DISTRIBUTION STUDY 
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Fig. 16 Transmitted noise levels into a double wail 
aluminum shell for different input locations 



Fig. 17 Natural frequencies of double wall composite shell 
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Fig. 18b Composite shell deflection levels for point 
loads acting on the exterior surface 
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Fig. 19 Composite shell deflection levels for point 
loads acting on the interior surface 
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Fig. 21 Sound pressure levels for aluminum and 
composite shells (exterior point loads) 




Fig. 22 Sound pressure levels for different structural and 
acoustic damping conditions (exterior point loads) 
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Fig. 23 Sound pressure levels in a composite shell 
for different loading conditions 




Fig. 24 Sound pressure levels in a composite shell 

for different fiber orientations (interior point 
loads ) 
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Fig. 27 Interior sound pressure levels for different shell 
configurations (no stringers) 
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Fig. 28 Interior sound pressure levels for different shell 
skin thicknesses 
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Fig. 29 Interior sound pressure levels for a shell with and 
without stringers 





Fig. 30 


Sound pressure levels for different locations of 
input forces 
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Sound pressure levels for a cylinder stiffened with 
heavy and light frames (ten frames) 
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Fig. 32 Interior sound pressure levels for a shell stiffened 
either with 10 heavy frames or 41 small frames 
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Fig. 33 Interior noise levels generated inside Skylab from 
operation of two different fans 
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Fig. 38 Sound pressure levels due to end plate 

vibrations with different loading conditions 
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Fig. 39 Sound pressure levels due to vibrating shell 
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Fig. 40 Total interior noise due to 120 dB exterior 
pressure 
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Fig. 41 Sound pressure levels due to point load inputs 
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Fig. 42 Total interior noise due to point load inputs 
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Fig. 50 Bending moment spectral densities for a cantilever box beam 
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Fig. 51 Deflection response spectral densities to calculated bending moment in 
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Fig. 53 Deflection response spectral densities for a coupled skin-stringer - beam 
and for a calculated bending moment input. 
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Fig. 54 Deflection response spectral densities for different geometric characteris- 
tics of the cantilever box beam. 
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Fig. 55 Deflection response spectral densities for the system shown in fig. 8 where 

(1) mi = k T i = k Rl =0 i 

(2) mi 7 ^ 0 , kn = k R i = 0 ; 

(3) k Ri ^ 0, mi = k Ti = 0 ; 

(4) k Ti ± 0, mi = k R i = 6 ; 

(5) mi ^ 0 , k T i ^ 0 , km ± 0 . 
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Fig. 60 Sound pressure levels in the interior for coupled and uncoupled conditions 
between the substructure and the noise transmitting stiffened sidewall. 
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Fig. 63 One-third octave sound pressure levels 
for a panel located at z = 0, x = 32 in 
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Sound pressure levels for vibrating panel 
located at z = 0, x = 10 in., y = 10 in. 
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Sound pressure levels for vibrating panel 
located at z = 0, x = 32 in., y = 10 in. 
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ia. 69 One-third octave sound pressure levels _ 

for a panel located at z = 20 in . , x - 84 m 
v = 198 in. 
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One-third octave sound pressure levels for a 
panel with dimensions L x = 21 in. and L 
= 26.1 in. y 
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Fig. 76 


One-third octave sound pressure levels for 
unstiffened panel 
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Fig. 78 


One— third octave sound pressure levels for 
a stiffened panel 





